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Abstract:  Optical  nanofibers  confine  light  to  subwavelength  scales, 

and  are  of  interest  for  the  design,  integration,  and  interconnection  of 
nanophotonic  devices.  Here  we  demonstrate  high  transmission  (>  97%) 
of  the  first  family  of  excited  modes  through  a  350  nm  radius  fiber,  by 
appropriate  choice  of  the  fiber  and  precise  control  of  the  taper  geometry. 
We  can  design  the  nanofibers  so  that  these  modes  propagate  with  most  of 
their  energy  outside  the  waist  region.  We  also  present  an  optical  setup  for 
selectively  launching  these  modes  with  less  than  1%  fundamental  mode 
contamination.  Our  experimental  results  are  in  good  agreement  with  simula¬ 
tions  of  the  propagation.  Multimode  optical  nanofibers  expand  the  photonic 
toolbox,  and  may  aid  in  the  realization  of  a  fully  integrated  nanoscale  de¬ 
vice  for  communication  science,  laser  science  or  other  sensing  applications. 
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1.  Introduction 

Optical  nanofibers  with  a  waist  diameter  smaller  than  the  wavelength  of  the  guided  light  are  cur¬ 
rently  used  for  non-linear  optics,  atomic  physics,  sensing,  and  fiber  coupling  [1-3].  The  intense 
evanescent  field  outside  of  an  optical  nanofiber  is  particularly  interesting  for  atom  trapping 
and  strong  atom-photon  coupling.  Tapering  standard  optical  fiber  down  to  submicron  diame¬ 
ters  has  been  a  successful  fabrication  technique  for  a  variety  of  applications  [4—8],  in  which 
transmissions  of  the  fundamental  mode  of  the  waveguide  reaching  more  than  99%  have  been 
achieved  [4,9].  Until  now,  the  efficient  guidance  of  higher-order  modes  has  not  been  observed, 
due  to  the  ease  with  which  they  couple  to  other  modes,  leading  to  large  losses.  This  restricts 
most  work  with  nanofibers  to  the  single  mode  regime,  where  the  diameter  is  small  enough  to 
only  support  the  fundamental  mode  HE\\  [3, 10], 

This  paper  reports  measurements  with  nanofibers  using  the  first  excited  TEq\,  TMq\,  and 
HE2 1  modes,  which  have  azimuthal,  radial,  and  hybrid  polarization  states,  respectively.  Prior 
work  on  higher-order  mode  propagation  has  shown  «30%  transmission  of  this  LP\  \  family  of 
modes  [11],  with  a  mode  purity  at  the  fiber  output  of  ss70%,  which  gives  a  total  transmission  of 
only  20%.  Here,  by  carefully  controlling  the  taper  geometry,  and  by  choosing  a  commercially 
available  50  micron  reduced-cladding-diameter  fiber  with  increased  numerical  aperture,  we 
demonstrate  transmissions  greater  than  97%  for  light  of  780  nm  wavelength  in  the  first  excited 
LP\  1  family  of  modes  through  fibers  with  a  350  nm  waist  radius.  Furthermore,  we  present  a 
setup  that  enables  us  to  efficiently  launch  these  three  modes  with  exceedingly  high  purity  at  the 
output,  where  less  than  1%  of  the  light  is  coupled  to  the  fundamental  mode  [12].  We  follow 
the  work  of  [13]  to  fabricate  the  nanofibers  and  use  a  series  of  diagnostics  during  the  pull  to 
monitor  the  quality  of  the  fiber.  We  record  the  total  transmission  of  the  light  and  image  the 
mode  exiting  the  fiber  for  the  duration  of  the  pull.  Analysis  of  the  transmission  as  a  function  of 
time  for  different  types  of  fibers  allows  us  to  estimate  which  modes  are  excited  during  the  pull 
as  well  as  their  relative  energies  of  excitation  through  the  use  of  spectrograms  [14]. 

These  results  open  the  way  to  efficiently  use  higher-order  modes  in  optical  nanofibers.  Unlike 
the  HE\  \  mode,  higher-order  modes  experience  a  cutoff  at  a  finite  radius.  This  allows  improved 
control  of  the  evanescent  field  extent  at  large  radii,  enabling  stronger  fibers  and  improved  han¬ 
dling  characteristics.  Our  work  enables  the  usage  of  higher-order  modes  for  atomic  physics 
applications.  In  particular,  the  spatial  interference  between  several  of  those  modes  can  create 
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unique  evanescent  field  distributions  on  the  waist,  providing  an  easy  and  self-consistent  way  to 
break  the  symmetry  along  the  propagation  axis,  suppressing  the  need  to  create  a  standing  wave. 
This  is  particularly  relevant  to  atomic  physics  applications  for  the  realization  of  a  one  color, 
blue-detuned  and  state  insensitive  trapping  potential  for  atoms  [15]. 

2.  Mode  propagation  in  an  optical  fiber 

References  [16, 17]  describe  the  modes  in  a  cylindrical  waveguide  using  Maxwell  equations. 
The  modal  fields  vary  as  exp  [Hfiz  —  Of)]  where  [3  is  the  propagation  constant  of  the  mode.  A 
mode  propagates  in  the  fiber  with  an  effective  index 


tleff  =  P/k.  (1) 

where  k  =  2n/\  is  the  free-space  propagation  constant  and  A  is  the  free-space  wavelength  of 
the  light.  The  propagation  of  light  inside  a  two-layered  step-index  fiber,  consisting  of  a  core  of 
radius  a  and  refractive  index  ncore  surrounded  by  a  cladding  of  radius  R  and  refractive  index 
ndad ,  depends  on  V, 


V  = 


2  na 


clad  ’ 


(2) 


V  plays  an  important  role  in  our  tapers,  since  the  radius,  a  in  Eq.  (2),  varies  enough  that  the 
interfaces  seen  by  the  modes  change  as  they  propagate  through  the  taper.  At  the  beginning  of 
the  taper,  the  light  is  confined  in  the  core,  and  guided  by  the  core-to-cladding  interface  with 
Vcore  as  in  Eq.  (2),  and  nciaij  <  neff  <  ncore.  At  the  end  of  the  taper,  the  core  is  negligible 
( cicore  ~10  nm  <C  A)  and  the  light  is  guided  by  the  cladding-to-air  interface  with  nalr  <  neff  < 
n clad-  Between  these  two  regimes,  the  light  escapes  the  core  to  the  cladding,  and  the  relevant 
radius  in  Eq.  (2)  becomes  R ,  which  is  much  larger  than  a.  Due  to  this  radius  increase,  and  the 
large  index  difference  between  and  na!>,  V  >>  1.  The  fiber  becomes  highly  multimode,  as 
the  number  of  bound  modes  is  proportional  to  V2/2  [16].  Maintaining  adiabaticity  through  this 
transition  is  critical. 


2.1.  First  family  of  excited  modes 

Figure  1(a)  shows  the  effective  index  of  refraction  against  V  for  several  low -order  modes  in  a 
nanofiber.  When  V  <  2.405,  the  fiber  supports  only  the  HE\\  mode.  For  a  typical  nanofiber, 
n a, re  ~  1.5  and  nciad  =  1.  The  TEq\  and  TMq\  modes  are  allowed  for  V  >  2.405,  and  the 
HE21  mode  is  allowed  for  V  >  2.8.  As  long  as  V  remains  lower  than  3.8,  only  the  four  modes 
mentioned  above  are  allowed.  In  the  weakly  guided  regime,  the  modes  of  interest  are  known 
as  the  LPn  family.  We  initially  launch  into  the  actual  LP\\  family,  because  at  the  entrance, 
the  fiber  is  weakly  guided.  The  taper  takes  us  into  the  strong  guiding  regime  where  the  LP\  \ 
splits  into  TE01,  TMqi.  and  II If  \  .  We  will  refer  to  the  set  of  modes  (77:: 01 .  TMq\.  and  HE21) 
during  the  entirety  of  the  pull  as  the  LP\  \  family.  This  simplifies  discussions  when  referring  to 
the  full  set  of  modes,  especially  in  reference  to  excitations  to  modes  or  families  with  the  same 
symmetry,  i.e.  LP\2  for  TE02,  TMq 1,  and  7/7si2-  Prior  work  has  emphasized  propagation  of  the 
HE  11  mode,  where  V  <  2.405  [3,7, 10].  We  are  interested  in  selectively  exciting  and  guiding 
the  LPn  family  through  a  nanofiber,  in  a  regime  where  2.405  <  V  <  3.8.  When  expanded  to 
free  space,  these  modes  have  the  intensity  and  polarization  profiles  shown  in  Fig.  1(b). 

2.2.  Adiabaticity  condition 

During  the  pull  we  continuously  decrease  R  and  the  light  escapes  the  core  to  the  cladding  as  it 
propagates  through  the  taper.  After  light  escapes  the  core,  the  presence  of  the  core,  the  cladding. 
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(a) 


TE„,  InO;  ™.,  HOC  HE ^O/- 


Fig.  1.  (a)  nefj  indices  of  several  low-order  modes  in  a  nanofiber  with  nciaci  =  1.5,  sur¬ 
rounded  by  vacuum  (nal>  =  1).  Below  V  ~  3.8,  two  families  of  modes  exist.  In  this  work, 
we  are  emphasizing  the  LP\  \  family  (circled).  By  symmetry,  modes  in  this  family  can  in¬ 
terfere  with  the  TMq2,  TEq2 ,  HE22  family,  which  may  be  excited  through  non-adiabatic 
processes,  (b)  Intensity  and  polarization  profiles  of  the  LP\\  family  of  modes  considered  in 
this  work. 


and  the  air  influence  the  mode.  The  fiber  is  highly  multimode,  and  modes  of  the  same  symmetry 
can  couple  to  each  other.  The  mode  evolution  in  a  taper  is  strongly  related  to  the  shape  of  the 
taper.  If  a  taper  is  too  steep,  the  mode  evolution  is  non-adiabatic  resulting  in  low  transmission. 
As  the  tapering  angle  Q.  is  reduced,  the  mode  propagation  becomes  more  adiabatic.  Following 
this  reasoning,  an  adiabaticity  criterion  has  been  derived  [16]  relating  the  characteristic  length 
of  the  taper  z.t,  to  the  characteristic  beat  length  between  two  modes  Zb  where  it  =  R/tan(f2), 
and  zb  =  2k /{Pi  -  fh)  =  A/(i neffl  -  neff2). 

The  mode  evolution  in  a  taper  is  adiabatic  when  the  fiber  is  long  enough  to  satisfy  Zt  >>  Zb 
[16].  For  a  typical  silica  optical  fiber  single  mode  at  A  =  780  nm,  we  calculate  using  this 
criterion  a  minimum  £2  of  a  few  milliradian  for  the  limiting  case  of  Zt  =  Zb-  This  implies  that 
a  taper  requires  sub-milliradian  Q.  to  achieve  adiabaticity  in  the  region  where  light  leaves  the 
fiber  core  and  becomes  a  cladding  mode.  Therefore,  choosing  an  optical  fiber  with  a  small 
initial  cladding  radius  (and  thus  a  reduced  Vciaci )  combined  with  a  high  numerical  aperture  is 
highly  advantageous  to  maintaining  adiabaticity  [11],  Additionally,  such  a  fiber  reduces  the 
overall  drawing  time,  length  requirements  of  the  pulling  apparatus,  and  the  overall  length  of  the 
taper. 

3.  Experimental  setup 

Figure  2  shows  a  diagram  of  our  experimental  setup.  Here,  we  review  the  mode  preparation, 
the  pulling  process,  and  the  detection  and  analysis.  We  efficiently  excite  the  TMq\ ,  TEq\,  and 
HEi\  modes  using  a  fiber-based  Cylindrical  Vector  Beam  (CVB)  generation  method  [12, 18]. 
The  fiber  is  drawn  using  a  heat-and-pull  method  [13, 19,20].  We  measure  the  transmission  while 
pulling,  recording  simultaneously  the  output  of  the  photodiode  to  a  digital  storage  oscilloscope 
and  beam  profiles  on  a  Charge-Coupled  Device  (CCD)  camera. 

A  New  Focus  Vortex  laser  delivers  a  TEMqq  Gaussian  beam  at  A  =7 8 0.24  nm.  This  wave¬ 
length  corresponds  to  the  D2  line  of  Rubidium,  our  atom  of  interest.  We  spatially  filter  this 
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Fig.  2.  (a)  Simplified  diagram  of  the  experimental  setup.  A  Gaussian  beam  passes  through 
a  ;r-phase  plate  and  is  coupled  into  the  fiber  to  be  drawn.  On  the  output  of  the  fiber,  a 
photodetector  (PD)  and  camera  (CCD)  monitor  the  transmission.  Typical  beam  images  are 
shown,  (b)  Schematic  of  tapered  nanofiber. 


beam  using  a  polarization-maintaining  optical  fiber,  and  collimate  with  an  asphere  to  a  1/e2  di¬ 
ameter  of  630  microns.  The  Gaussian  beam  passes  through  a  phase  plate  that  imparts  a  n  phase 
shift  on  half  of  the  beam,  producing  a  two-lobed  beam  that  approximates  a  TEMq\  free-space 
optical  mode.  The  phase  plate  is  on  a  translation  stage,  so  that  we  can  easily  switch  between  the 
LP\\  and  LPoi  families  without  adjusting  alignment.  Figure  2  shows  profiles  of  the  beam  be¬ 
fore  and  after  the  phase  plate.  We  couple  the  beam  into  the  fiber  using  a  matched  asphere.  The 
coupling  coefficients  to  the  fiber  modes  are  determined  by  the  beam  symmetry.  The  inversion 
of  the  polarization  over  half  the  incident  beam  allows  us  to  selectively  excite  the  LP\  i  family 
with  high  efficiency.  Stress-induced  birefringence  at  the  input  end  of  the  fiber  can  be  used  to 
selectvely  excite  the  TMq\ ,  TEq\,  and  ///A]  modes  individually  [12],  but  we  have  not  used  that 
selectivity  here.  The  fundamental  HE\  |  mode  is  only  excited  through  aberrations  in  the  beam. 
To  achieve  efficient  coupling  into  the  LP\  \  family,  we  wrap  the  fiber  with  two  or  three  wind¬ 
ings  around  a  4-mm  diameter  mandrel  that  attenuates  any  modes  higher  than  the  LP\  i  family. 
Figure  2(a)  shows  typical  input  (black  and  white)  and  output  (color)  modes. 

We  detail  our  nano  fiber  pulling  process  in  Ref.  [14,20].  We  start  by  clamping  a  fiber  to  two 
high-precision  Newport  XML  210  motor  stages.  We  strip  the  buffer  from  the  fiber  and  clean 
it  thoroughly  following  the  procedures  in  [20].  We  image  the  fiber  with  an  optical  microscope 
and  remove  any  undesired  particle  remaining  on  the  fiber  surface.  This  is  important  for  the 
repeatability  of  our  measurements.  An  oxyhydrogen  flame  with  a  stoichiometric  combination 
brings  a  0.75  mm  long  portion  of  the  fiber  to  a  temperature  that  exceeds  the  softening  point  of 
fused  silica.  We  pull  on  the  fiber  ends  at  a  typical  relative  velocity  of  0.1  mm/s.  We  calculate 
the  stage  trajectories  to  produce  a  fiber  of  a  chosen  geometry  using  an  algorithm  that  relies 
on  conservation  of  volume  [5].  We  divide  the  pull  into  approximately  100  steps.  Each  step 
adds  a  small  section  to  the  taper  that  reduces  the  radius  of  the  waist.  The  sections  are  small 
enough  to  be  considered  linear  locally.  The  compilation  of  each  small  taper  creates  the  final 
taper  with  a  desired  geometry,  which  is  generally  composed  of  a  few  mrad  steep  linear  section 
that  reduces  to  a  radius  of  6  /xm,  and  then  connects  to  an  exponential  section  that  gradually 
reaches  a  submicrometer  radius.  The  central  region  forms  a  7-mm  uniform  waist  (see  Fig.  2(b) 
for  the  geometry).  To  observe  all  possible  mode  cutoffs,  the  fibers  are  typically  drawn  to  R  «280 
nm,  at  which  point  only  the  fundamental  ///-)  |  mode  propagates.  We  vary  £1  from  0.4  mrad  to  4 
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mrad,  resulting  in  pull  times  lasting  between  100  to  1000  seconds  The  output  side  of  the  fiber  is 
held  straight  with  no  mode  filtering.  We  follow  the  mode  evolution  for  the  duration  of  the  pull 
by  monitoring  the  transmission  of  a  few  mW  of  laser  power  through  the  fiber.  The  transmitted 
beam  is  monitored  both  by  a  CCD  and  by  a  photodetector  (PD).  The  beamsplitter  shown  in  the 
figure  is  tilted  to  as  small  an  angle  as  possible  to  eliminate  polarization-dependent  reflections. 
Another  PD,  not  shown  in  Fig.  2(a),  measures  the  input  laser  power  during  the  pull  to  normalize 
the  transmission  signal.  We  record  PD  signals  with  a  Tektronix  DPO7054  oscilloscope  with  16- 
bit  resolution  and  a  sample  rate  of  1-10  ksamples/s. 

4.  Results 

In  this  section,  we  present  our  results  on  the  transmission  of  the  LPi  i  mode  family  in  the  context 
of  adiabaticity.  We  first  analyze  the  improvements  gained  by  choosing  a  fiber  with  a  reduced 
cladding  radius  and  increased  numerical  aperture.  In  [11],  Frawley  et  al.  looked  at  the  improve¬ 
ment  obtained  by  reducing  the  initial  R  from  125  pm  to  80  /tin.  Moving  from  80  p m  to  50  p m 
fibers,  we  improve  the  transmission  from  10%  to  51%  for  a  2-mrad  taper.  A  spectrogram  anal¬ 
ysis  depicts  fewer  and  weaker  excitations  of  higher-order  modes.  Second,  we  demonstrate  an 
improvement  in  adiabaticity  through  the  control  of  the  taper  geometry.  Our  ability  to  vary  £2 
allows  for  marked  improvement  in  guidance  efficiency  through  the  nanofiber,  increasing  from 
16%  at  4  mrad  up  to  97.8%  at  0.4  mrad. 

4.1.  Varying  the  fiber  type 

4.1.1.  Transmission  measurements 

Figure  3  shows  typical  transmissions  obtained  when  exciting  the  fiber  with  a  combination  of 
LP\  |  modes.  We  plot  the  normalized  transmission  during  a  single  pull  as  a  function  of  the  fiber 
waist  radius.  We  can  identify  four  distinct  regimes  in  Fig.  3.  The  modes  are  initially  confined  to 
the  core  (regime  A  in  Fig.  3).  Adiabaticity  can  easily  be  achieved,  and  the  transmission  is  steady. 
Because  this  fiber  initially  also  supports  the  LP02  and  Ll\\  modes  and  because  our  launch  can 
weakly  excite  these  modes,  we  observe  a  slight  drop  in  transmitted  power  near  R  =  20  pm, 
when  these  higher-order  modes  become  cladding  modes.  Regime  B  occurs  after  the  light  has 
escaped  from  the  core  to  the  cladding.  Because  the  cladding  is  typically  much  larger  than  the 
core  (R/a  >  10),  and  because  nciad  —  najr  ncore  —  nciad,  V  increases  by  over  two  orders 
of  magnitude  (Vciad  ~  200).  If  the  core-cladding  transition  is  not  adiabatic,  modes  of  similar 
symmetry  are  excited.  In  particular,  we  observe  transfers  of  energy  to  the  LP\2  family  which 
contains  the  modes  TEq2,  'I'M 02,  and  lllfii-  The  interaction  between  those  modes  results  in 
mode  beating  inside  the  fiber,  and  an  oscillation  in  the  amount  of  output  light.  The  oscillations 
continue  to  R  ss  0.7  pm  (regime  C),  where  for  a  typical  silica  fiber,  Vciatj  ss  6.  From  this  cutoff 
location,  it  is  clear  that  much  of  the  beating  behavior  is  due  to  the  LP12  modes  which  have 
been  excited  through  non-adiabatic  transitions.  The  pull  extends  through  the  TMq\,  TEq\ ,  HE21 
cutoffs  (regime  D ).  After  reaching  the  cutoff  radius  near  R  =  290  nm,  very  little  light  reaches 
the  photodetector  -  typically  less  than  1-2%  -  indicating  low  population  of  the  fundamental 
mode.  To  determine  the  loss  in  the  HE\\  mode,  we  measured  the  transmitted  power  before 
and  after  the  pull  by  translating  the  phase  step  out  of  the  way.  The  losses  in  this  mode  were 
negligible,  meaning  that  we  reach  a  mode  purity  of  98-99%.  Mode  purity  is  of  fundamental 
importance  for  nanofiber  applications,  and  in  particular  for  atomic  physics  where  polarization 
control  and  stability  are  required. 

We  have  compared  the  transmissions  obtained  using  fibers  of  diameter  80  pm  and  50  pm, 
and  observed  the  beneficial  effects  of  smaller  clad  fibers  and  higher  numerical  aperture  for 
£2  =  2  mrad.  This  taper  angle  is  chosen  because  it  is  non-adiabatic  and  highlights  the  effect 
of  fiber  diameter  on  adiabaticity.  The  improvement  in  adiabaticity  is  clear  in  Fig.  3,  with  a 
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Fig.  3.  Evolution  while  tapering  of  the  transmission  through  fibers  with  a  half  angle  of  2 
mrad  as  a  function  of  the  radius  of  the  waist,  (a)  Fiber  with  an  initial  diameter  of  80  pm.  (b) 
Fiber  with  an  initial  diameter  of  50  pm.  The  spectrograms  associated  with  those  transmis¬ 
sion  curves  give  a  clear  picture  of  the  power  transfers  during  the  pull.  Note  the  logarithmic 
scales  on  the  horizontal  axis,  (c)  (1.0  MB)  Movie  of  the  evolution  of  the  beam  transmitted 
through  the  fiber  measured  on  the  CCD  during  a  2  mrad  pull  of  the  50  p m  fiber.  We  record 
one  frame  every  second,  and  display  them  at  a  7  frames  per  second  speed.  Sections  A,  B, 
C,  and  D  are  described  in  the  text.  The  properties  of  the  fibers  used  are  summarized  in  the 
table. 


transmission  of  52%  obtained  for  the  50-/im  fiber,  compared  to  the  10%  transmission  of  the  80- 
jum  fiber.  The  table  in  Fig.  3  compares  the  properties  of  the  80- pm  and  50-/im  fibers.  Note  that 
the  SM1500  fiber  initially  supports  the  LPq\,  LP\\  as  well  as  the  next  families  of  excited  modes 
LPq2  and  LPi  \  (VCOre,R0  =  4.3).  These  higher-order  modes  are  substantially  filtered,  though  not 
completely,  by  winding  the  fiber  around  a  4  mm-diameter  rod. 

Using  the  numerical  aperture  (NA)  and  the  cutoff  wavelength  provided  by  Thorlabs  and 
Fibercore,  we  are  able  to  derive  using  Eq.  (2)  the  initial  V  value  Vcore,R0  ■  the  radius  at  which  the 
LP\  \  family  escapes  from  the  core  to  the  cladding  R1/’1 1  (VcoreRc  =  2.405:  the  mode  enters  the 
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cladding-to-air  guidance  regime)  and  the  corresponding  fiber  V -number  at  that  radius  Vci,uI.rc 
for  a  wavelength  of  780  nm.  Empirically,  we  obtain  an  estimate  of  R^Pn  on  the  transmission 
plot  when  the  oscillations  start.  The  numbers  agree  quantitatively.  We  see  that  reducing  the 
initial  cladding  radius  and  increasing  the  fiber  numerical  aperture  allows  the  LP\\  family  to  leak 
from  the  core  to  the  cladding  at  a  smaller  fiber  radius.  This  results  in  a  significantly  reduced 
Vdad  when  the  modes  escape  from  the  core  to  the  cladding:  By  reducing  the  cladding  diameter 
from  80  ji  m  to  50  jum  and  increasing  the  numerical  aperture  from  0.18  to  0.30,  the  number 
of  available  modes  decreases  by  more  than  an  order  of  magnitude  for  these  two  commercially 
available  fibers.  Improvements  could  be  achieved  by  pre-etching  the  fiber  to  a  smaller  diameter 
so  that  the  initial  core  radius  to  cladding  radius  ratio  R/a  is  further  reduced.  In  this  case,  the 
numerical  aperture  remains  unchanged,  and  the  number  of  available  modes  directly  scales  with 
the  square  of  the  ratio  of  initial  cladding  radii. 

4.1.2.  Spectrograms 

In  non-adiabatic  propagation,  the  LP\  \  modes  couple  to  higher-order  modes  of  the  same  sym¬ 
metry,  belonging  to  families  LP\m  (m  >  2).  Because  they  propagate  with  different  propagation 
constants  during  the  pull,  they  accumulate  a  phase  difference  leading  to  interference  in  the 
amount  of  light  that  coupled  back  into  the  core.  Since  the  photodetector  only  measures  core 
light,  this  interference  leads  to  oscillations  in  the  transmission  (Fig.  3).  A  useful  way  to  exam¬ 
ine  these  data  is  through  spectrograms,  which  plot  local,  windowed  Fourier  transforms  of  the 
transmission  signals  as  a  function  of  waist  radius.  Plotting  the  spectrogram  of  the  transmission 
signal  [14,21],  we  directly  observe  the  contribution  of  various  pairs  of  modes  in  the  beating. 
Figure  3  shows  the  spectrograms  for  the  80-/2 m  and  50-/1  m  diameter  fiber  pulls.  Each  line  in 
the  spectrogram  is  specific  to  the  beating  between  a  LPoi  or  LP\  \  mode  and  another  mode  of 
similar  symmetry  excited  during  the  pull.  The  curve  ends  when  one  of  the  modes  reaches  its 
cutoff:  the  energy  is  then  lost  via  coupling  to  radiative  modes. 

The  number  of  lines  observed  in  a  spectrogram  is  directly  related  to  the  excitation  of  higher- 
order  modes  through  non-adiabatic  processes  from  a  single  launched  mode:  the  more  lines 
present  in  the  spectrogram,  the  less  adiabatic  the  pull  is.  Moreover,  the  colormaps  in  Fig.  3  are 
normalized  in  such  a  way  that  the  intensity  of  each  red  line  gives  the  strength  of  the  energy 
transferred.  It  is  clear  from  Fig.  3  that  more  intense  lines  are  present  in  the  80-/2 m  fiber  than  in 
the  5 0-/rm,  further  supporting  our  observation  of  more  stringent  adiabaticity  requirements  for 
fibers  with  a  large  cladding  radius  and  small  numerical  aperture. 

4.1.3.  Imaging  the  fiber  output 

We  monitor  the  transmitted  beam  with  a  CCD,  and  obtain  Media  1  showing  both  the  core 
and  cladding  light  as  a  function  of  time  throughout  the  tapering  process.  Using  a  microscope 
objective,  we  first  image  the  end  of  the  fiber  to  observe  the  core-guided  light  (Media  1  in  Fig.  3). 
The  movie  shows  oscillations  that  result  from  mode  competition,  which  modulates  the  amount 
of  light  that  exits  the  fiber  in  the  core. 

Media  2  in  Fig.  4(a)  also  depicts  the  transfer  of  energy  from  the  core  to  the  cladding.  Because 
the  cladding  intensity  is  low  compared  with  the  core,  we  record  two  spatially  separated  images 
of  the  fiber  simultaneously  by  using  strong  and  weak  reflections  from  a  thick  beamsplitter  with 
AR-coated  front  face.  Figure  4  shows  the  normalized  fractions  of  energy  exiting  the  fiber  from 
the  cladding  and  from  the  core.  The  two  signals  are  out  of  phase.  We  note  that  the  sum  of 
the  energy  contained  in  the  core  and  the  energy  contained  in  the  cladding  does  not  add  up 
to  the  total  energy  input  in  the  fiber.  Outside  the  taper  region,  cladding  light  becomes  highly 
scrambled  and  lost  through  the  fiber  buffer,  resulting  in  both  the  speckle  observed  in  Media  2 
and  the  reduced  total  transmitted  power.  Residual  cladding  light  is  spatially  filtered  from  hitting 
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Fig.  4.  Amount  of  light  (normalized)  exiting  the  fiber  from  the  core  (blue  curve)  and  from 
the  cladding  (red  curve).  The  signals  are  out  of  phase,  confirming  the  transfer  of  energy 
between  modes  during  the  tapering.  We  observed  the  two  simultaneously  by  using  the  two 
reflections  from  a  thick  beamsplitter,  (a)-(b)  (1.1  MB)  Media  2  shows  the  evolution  of  the 
beam  transmitted  through  a  nanofiber  during  a  portion  of  a  pull,  where  the  power  is  high 
enough  to  observe  the  cladding  light. 


the  photodetector,  so  that  the  observed  oscillations  Fig.  3  are  due  only  to  core-guided  light. 

4.2.  Varying  the  angle 

The  measurements  in  this  section  use  the  reduced-cladding  Fibercore  SM1500.  For  this  fiber, 
the  LP\  i  modes  transition  to  cladding  modes  near  R  =  1 3  fim.  By  using  a  reduced-diameter 
fiber,  we  observed  a  drastic  improvement  of  the  transmission  of  the  LP\\  modes.  To  further 
improve  the  adiabaticity,  it  is  necessary  to  look  into  more  details  of  the  tapering  process  itself. 
We  have  lowered  £2  to  improve  the  transmission  over  what  is  observed  in  the  previous  section. 
Figure  5  shows  the  results  of  draws  using  £2  =  4,  2,  1,  0.75  and  0.4  mrad.  Although  each  plot 
shows  the  same  qualitative  behavior  as  described  in  Fig.  3  for  the  2-mrad  pull,  the  strength  of 
the  features  depends  on  £2. 

The  free-space  mode  at  the  fiber  input  has  a  spatial  polarization  that  is  an  equal  superposition 
of  HE2\  and  TMq\  (or  TEq\).  However,  mode  conversion  occurs  where  the  fiber  is  wound 
around  the  mode-filtering  mandrel  so  that  the  distribution  entering  the  nanofiber  is  unknown. 
Within  the  nanofiber  waist,  which  is  held  fixed  and  straight,  mode  conversion  is  unlikely  to 
occur  so  that  the  desired  LP\  \  mode  can  be  achieved  after  the  pull  [12].  By  R  =  0.45  /im,  only 
the  TMqi,  TEcn,  and  HE2\  modes  are  confined,  with  a  small  contribution  in  the  fundamental 
HEn.  The  II P\2  mode  achieves  cutoff  at  Vcja([  =  2.8  (R  ~  330  nm),  earlier  than  the  TMq\ 
and  TEq\  modes,  which  reach  cutoff  at  V  =  2.4  (R  «  290  nm).  At  R  sa  330  nm,  the  power  is 
reduced  by  the  HE2\  content,  which  is  determined  by  the  initial  superposition  of  states  entering 
the  nanofiber. 

The  transmitted  power  in  the  LP\\  mode  family  is  16.6%  for  £2  =  4  mrad.  The  transmitted 
power  drops  sharply  near  R  =  13  fim,  and  undergoes  strong  oscillations  between  10-80%.  For 
£2  =  2  mrad,  the  oscillations  below  /?  =  13  jX  m  are  reduced,  with  the  transmission  fluctuating 
between  30-90%,  and  the  transmitted  power  improves  to  5 1 .5%.  Further  decreasing  the  angle  to 
£2=1,  0.75  and  0.4  mrad,  improves  the  transmission  to  84.3%,  91.2%  and  97.8%  respectively, 
where  uncertainty  is  dominated  by  systematic  effects  that  should  be  less  than  1%.  Figure  5  also 
shows  excellent  agreement  between  the  experimental  results  and  those  obtained  using  commer¬ 
cial  waveguide  propagation  software  [22].  For  those  calculations,  the  propagation  was  modeled 
using  the  targeted  fiber  geometry.  In  our  transmission  measurements,  the  scaling  between  time 
and  radius  is  made  using  a  separate  algorithm  that  models  the  dynamics  of  the  pull  based  on 
conservation  of  volume  [20].  For  £2  =  0.4  to  2  mrad,  the  mode  cutoff  positions  we  obtain  using 
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Fig.  5.  (a)  Transmission  through  a  SMI 500  nanofiber  for  SI  =  4  mrad  (green),  2  mrad 
(blue),  1  mrad  (purple),  0.75  mrad  (yellow),  0.4  mrad  (red).  The  transmission  plots  have 
the  fundamental  mode  subtracted  out  (typically  about  1%)  to  accurately  describe  the  total 
transmission  of  the  LP\  i  modes.  The  horizontal  axis  for  the  4  mrad  pull  was  renormalized  to 
take  into  account  extra  tension  in  the  fiber  due  to  the  rapidity  of  the  pull,  (b)  Simulated  (red 
circles)  and  experimental  (blue  lines)  final  transmissions  through  the  fiber  as  a  function 
of  angle.  Decreasing  SI  enables  us  to  improve  the  transmission  of  the  LP\\  family  up  to 


97.8%. 


the  scaling  from  the  algorithm  correspond  directly  to  what  is  expected  theoretically.  We  also 
observe  excellent  agreement  between  the  experimental  transmission  measurements  and  the  cal¬ 
culations,  confirming  the  accuracy  of  our  pulling  procedure.  We  were  not  limited  by  systematic 
effects  that  might  include  accumulation  of  contaminants  for  longer  pulls,  and  asymmetric  pro¬ 
files  for  faster  ones.  For  the  fastest  pull  (4  mrad),  we  observed  a  discrepancy  due  to  systematic 
effects,  and  we  must  apply  a  different  scaling  to  match  the  theoretical  cutoffs. 

For  SI  =  0.4  mrad,  the  most  shallow  angle  studied,  the  amplitude  of  the  oscillations,  which  is 
directly  related  to  the  energy  transfer  to  undesired  modes,  is  reduced  to  a  few  percent.  The  ob¬ 
served  transmissions  are  due  to  losses  into  and  out  of  the  waist.  Because  the  fiber  is  symmetric, 
the  normalized  transmitted  power  is  the  square  of  that  in  the  waist.  For  SI  =  0.4  mrad,  this  leads 
to  98.9%  power  in  the  waist.  We  believe  that  such  a  fiber  is  usable  for  various  applications  in¬ 
volving  higher-order  modes.  We  note  that  reaching  adiabaticity  for  higher-order  modes  requires 
fibers  that  are  substantially  longer  than  for  the  fundamental  mode.  For  the  H E\  \  mode,  trans¬ 
missions  greater  than  99%  can  easily  be  achieved  for  SI  up  to  5  mrad  [14].  We  also  observed 
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that  when  we  remove  the  /r-phase  plate  from  the  launch,  the  transmission  in  the  fundamental 
mode  is  essentially  equal  to  the  transmission  before  pulling,  as  adiabaticity  is  strongly  satisfied 
for  this  mode.  During  the  tapering  process,  the  higher-order  modes  escape  from  the  core  to 
the  cladding  earlier  than  the  fundamental  HE\  \  mode.  When  the  HE\  \  mode  finally  transitions 
to  a  cladding  mode,  R  has  decreased,  so  that  V  and  the  number  of  available  modes  to  couple 
to  is  smaller.  The  reduction  in  R  also  leads  to  an  increase  in  the  difference  between  adjacent 
propagation  constants,  allowing  less  mode  interaction  and  a  steeper  £2.  Mode  conversion  also 
occurs  throughout  the  fiber  and  not  just  at  the  core-cladding  transition  point,  but  this  region  has 
the  most  stringent  adiabatic  criterion. 

5.  Conclusion 

We  have  demonstrated  propagation  of  higher-order  modes  in  nanofibers  using  the  TEq\ ,  TMq\, 
and  HEj \  modes.  By  tapering  the  fiber  with  angles  near  0.4  mrad  and  using  a  commercial,  off- 
the-shelf  fiber  with  50  /um  diameter,  we  have  achieved  transmission  efficiency  of  97.8%  with 
excellent  mode  purity,  a  factor  of  four  higher  than  previous  work,  and  more  than  one  order  of 
magnitude  improvement  on  mode  purity.  Critical  to  this  work  was  a  spectrogram  analysis  of 
the  modes  present  during  the  pulling.  Our  experimental  results  agree  with  simulations  of  the 
propagation  through  the  taper.  High  transmissions  of  LP\  \  modes  with  high  purity  is  a  promising 
tool  for  atomic  physics,  expanding  the  possible  intensity  and  polarization  configurations  of 
evanescent  fields  surrounding  the  nanofiber. 
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